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It is numerically investigated whether a deformable object can propel itself in a superfluid. Articulated bodies and 
multi-component condensates are examined as swimmers. An articulated two-body swimmer cannot obtain locomotion 
without emitting excitations. More flexible swimmers can do so without the need to excite waves. 


1. Introduction 

Superfluids have no viscosity. When an object is moved in 
a superfluid at a subcritical constant velocity, it experiences 
no drag force. It would seem, therefore, that the arms and legs 
of a swimmer will be unable to generate propulsion in a su¬ 
perfluid. Can we swim in such a fluid? In other words, is the 
self-propulsion of a deformable object in a superfluid possi¬ 
ble? 

The present paper shows that the answer to this question is 
yes. Let us first consider a case in which the velocity of the 
flow around the swimmer is much slower than the velocity 
of sound. In this case, a superfluid may be regarded as an in¬ 
compressible ideal fluid. Several authors have shown that the 
self-propulsion of a deformable body is possible even in an in¬ 
compressible ideal fluid, and therefore, we expect that self¬ 
propulsion is also possible in a superfluid. When the swimmer 
deforms quickly, on the other hand, the superfluid cannot be 
regarded as incompressible, and waves and quantized vortices 
are generated. Since these excitations are able to carry mo¬ 
mentum, the swimmer can obtain thrust by releasing them. 
The aim of the present paper is to verify these speculations. 

The dynamics of objects immersed in superfluids have been 
studied in various situations, such as using particles to visu¬ 
alize superfluid flow,^’^ using oscillating spheres^ and wires^ 
to disturb superfluid helium, and considering the behavior of 
an ion^^ in a gaseous Bose-Einstein condensate (BEC). How¬ 
ever, the self-propulsion of an object in a superfluid has not 
been studied. 

In the present paper, we use the mean-field Gross-Pitaevskii 
(GP) theory to numerically demonstrate that deformable ob¬ 
jects can swim in a two-dimensional (2D) superfluid. Two 
kinds of swimmers are examined: articulated bodies with 
movable joints and multicomponent BECs. Eor the latter type, 
the swimmer is also a BEC and can be deformed by chang¬ 
ing the interactions between atoms. Since a superfluid has no 
viscosity, one may think that the swimmer would behave as 
if it had a high Reynolds number. However, the condition for 
self-propulsion in a superfluid is shown to be similar to that 
for a swimmer with a low Reynolds number in a viscous fluid; 


this is reminiscent of Purcell’s “scallop theorem”. 

This paper is organized as follows. Section 2 investigates 
the dynamics of articulated bodies in a BEC. Section 3 studies 
the dynamics of multicomponent BECs. Section 4 presents 
the conclusions of our study. 

2. Articulated bodies in a superfluid 

We consider the dynamics of a superfluid of atoms with 
mass m in the mean-field approximation. The dynamics of the 
macroscopic wave function are described by the GP equa¬ 
tion 

= -LvV + vip + (1) 

ot 2m 

where V is the potential produced by a swimmer. Eor sim¬ 
plicity, we will consider a 2D space, which can be realized 
by tightly confining the system in the z direction. Assuming 
the form of the wave function is 0 ^ 3 D(x,y, z) = 0^(x,y)0(z), 
the effective interaction coefficient g in Eq. (1) is given by 
g = 4nffa f \(l)\^dzlm, where a is the v-wave scattering length. 

In this section, we consider swimmers that consist of el¬ 
liptic bodies, as illustrated in Eig. 1. The elliptic bodies are 
connected at joints located on the long axes of the ellipses. 
The distance between the centers of the ellipses and the joints 
is d. The angles x of the joints are controlled variables. To 
reduce numerical errors arising from the rigid boundary con¬ 
dition at the surface of the bodies, a “soft” boundary condition 
is adopted. Eor the jth elliptic body centered at (xj, yj) and for 
which the angle between its long axis and the x axis is 6j, the 
potential is given by 

Vj(x,y) = HQ exp{/pot[l - {^jlaf - {rjjHif]], (2) 

where = {x-xf) cos 6j -f iy-yj) sin 6j, rjj = -(x-xj) sin 6j + 
(y - yj) cos 6j, and po is the chemical potential far from the 
swimmer. The potential in Eq. (2) is vj = po at an ellipse 
-f (jijlP)^ = 1, which corresponds to the surface of the 
body. The parameter /pot in Eq. (2) determines the sharpness 
of the surface; the Gaussian potential rises steeply for a large 
/pot- The results are insensitive to the value of /pot, and we set 
/pot = 1 in the following calculations. The potential V in the 
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GP equation is thus given by 

y{x,y) = ^Vj{x,y). (3) 

j 

The Lagrangian for the swimmer is given by 

^ = Z [y^4) 


where M and I are the mass and the moment of inertia of 
the elliptic body, respectively. The interaction potential f/pot 
between the swimmer and superfluid has the form 

Up„,{X,Y,9,x\^) = f V\4>Ur, (5) 

where (X, Y) are coordinates of the center of mass (COM) of 
the swimmer. The potential V depends on (X, T), 6, and 
through Eq. (2). In the following, we assume M = najSmno 
and I = M{a^ + 0^)14, where no is the density far from the 
swimmer. 

Equation (1) and the equation of motion for the swimmer 
are solved numerically, using the pseudo spectral and fourth- 
order Runge-Kutta methods. The initial state of j/r is pre¬ 
pared by the imaginary-time propagation method, in which 
the swimmer is at rest in the initial form. The density 
far from the swimmer is constant no, which gives the char¬ 
acteristic length ^ = filimgno)^^^ and time r = ^/v^, where 

= (g^olmy^^ is the velocity of sound. The following quan¬ 
tities are fixed: a = 5^, jS = and d = 1^. A periodic 
boundary condition is imposed by the pseudo spectral method. 
The numerical space is taken to be sufficiently wide, and 
the boundary does not affect the results if the motion of the 
swimmer is moderate. In this case, the total momentum of 
the swimmer and surrounding fluid is conserved. When the 
swimmer moves violently, waves and vortices are released, 
and these reach the boundary. In that case, an artificial damp¬ 
ing term -T(r)fidi//ldt is added to the left-hand side of Eq. (1), 
where r(r) is nonzero only near the boundary. The artificial 
damping term can suppress the disturbances at the boundary 
without affecting the dynamics of the swimmer. 


2 .1 Two-body swimmer 

We begin by considering a swimmer created from two el¬ 
liptic bodies connected by a joint, as shown in Eig. 1(a). The 
state of the swimmer is specified by the COM coordinate 
(X, T), direction 6 of the symmetry axis, and their time deriva¬ 
tives (X, Y), 6. The angleis the controlled variable, and the 
shape of the swimmer has only one degree of freedom. 

Eor the two-body swimmer, the Lagrangian in Eq. (4) is 
rewritten as 

L = M(X2 + Y^)+^ _ (^2 2xi] 

+1{0^ + - Upot(X, Y, 0,xi; <1^). (6) 

The dynamics of the swimmer are thus described by the Euler- 




Fig. 1. Schematic illustration of the two-body and three-body swimmers 
used in Secs. 2.1 and 2.2, respectively. The distance d is constant between the 
center of the elliptic bodies and the joints. The angles and;^2 are controlled 
parameters. 
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Fig. 2. (Color online) Time evolution of (a) the density profile and (b) 
the position X of the COM. The angle is driven as in Eq. (8) with A = 
7r/4 and cu = 0.025T“k The length and time are normalized by the healing 
length ^ and r = ^/v^, where is the sound velocity. See the Supplemental 
Material for a movie of the dynamics shown in (a).^^ 


Lagrange equations, 

2MX 

2MY 

[Md\l-cos2xi) + 2l]e 


dU, 


pot 


dx ’ 


dU, 


pot 


dY ’ 


(7a) 

(7b) 


dU, 


pot 


de 


-2Md Oxx sin2qfi. 

(7c) 


Equations (1) and (7) determine the dynamics of the system. 
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Fig. 3. (Color online) Time evolution of (a) the density profile |j/rp, and 
(b) the position X of the COM for oj = 0.4t“^ The angle is driven as 
in Eq. (8) for ojt < 20n, md xi(t) = 0 for cot > 20n. The dashed line in 
(b) indicates cot = 20n. See the Supplemental Material for a movie of the 
dynamics shown in (a).^^ 


First, we show that such an object cannot generate locomo¬ 
tion when it is deformed slowly. The swimmer is deformed 
periodically as 

/ 1 n (jjt\ 

where A12 is the amplitude and u; is a frequency of the os¬ 
cillation. Figure 2 shows the time evolution of the system for 
A = 711A and oj = 0.025t“^ The swimmer is initially formed 
by two elliptic bodies that are joined at an angle 7r/2, and then 
it closes and opens periodically, as shown in Fig. 2(a). We can 
see in Fig. 2(b) that the action of the swimmer generates no 
net locomotion, whereas the position X of the COM moves 
back and forth at a frequency oj, while conserving the total 
momentum. 

The fact that the swimmer cannot propel itself in Fig. 2 
reminds us of Purcell’s “scallop theorem”,^^ which states that 
when the motion of a swimmer has time-reversal symmetry, 
there is no net locomotion of the swimmer in a viscous Stokes 
fluid. The periodic motion in Eq. ( 8 ) has time-reversal sym¬ 
metry, and therefore the swimmer cannot swim in a Stokes 
fluid. Reference 15 showed that the scallop theorem can be 
generalized to a perfect fluid, i.e., a periodic motion with time 
reversal symmetry never generates locomotion in a perfect 
fluid. When the motion of the swimmer is sufficiently slow, 
and no waves or quantized vortices are excited, a superfluid 
can be regarded as a perfect fluid. The generalized scallop the¬ 
orem is therefore applicable to the situation shown in Fig. 2. 
Thus, self-propulsion in a superfluid is impossible for a two- 
body swimmer with sufficiently slow deformation. 

In the dynamics shown in Fig. 2, the typical velocity of 


the elliptic bodies is ojd = 0.175v^, and this is much smaller 
than the velocity of sound v^; thus, the fluid around the swim¬ 
mer is barely disturbed. When the swimmer is deformed more 
quickly, on the other hand, waves and quantized vortices are 
created, and the condensate can no longer be regarded as a 
perfect fluid; the generalized scallop theorem is no longer 
applicable. Figure 3 shows the dynamics of the system for 
OJ = 0.4t“\ which is 16 times larger than the value used for 
Fig. 2. In this case, the shape of the swimmer changes faster 
than the sound velocity, and waves and solitonic pulses^^’^^ 
are emitted from the swimmer, as shown in Fig. 3(a). By 
disturbing the fluid, the swimmer gains net locomotion, as 
shown in Fig. 3(b). At ojt = 20;: [the vertical dashed line in 
Fig. 3(b)], i.e., after ten strokes, the deformation of the swim¬ 
mer is stopped: Ad (0 = 0 ojt > 2071. Even after that, the 
swimmer continues to travel, which indicates that the swim¬ 
mer acquires momentum by emitting waves and solitonic ex¬ 
citations into the surrounding fluid. 


2.2 Three-body swimmer 

We next study the dynamics of a swimmer consisting of 
three elliptic bodies connected by two joints, as illustrated in 
Eig. 1(b). The Lagrangian in Eq. (4) can be written as 

Md^ 


3M .o .o 

- {X^ + Y^) + 

2 3 


+Ti)^ + 0 +T2)^ 


+(^ +Ti)(^ +T 2 ) cos(ai -X 2 ) + ^0(6 -i-Afi) cosAd 
-\-30(0 d-xi) COSX 2 + 


(9) 


’9 0 +Ti)^ + 0 +T 2 )^] 


-f/pot(2f, T,6>, An, Ar2;<A). 
which gives the Euler-Lagrange equations, 
dU^ot 


3MX 

3MY 


dX ’ 

dU, 


pot 


dY ’ 

5 t/pot Mcf- 


(10a) 

(10b) 

[(xi +X 2 ) (2 + cos(yi -Xi)) 


60 3 

-(20+xi +X 2 )(X\ -A" 2 )sin(yi -xi) 

+3x\ cos;(fi + 3x2 cosxi - 3(20 +xi)Xi sin;(fi 

-3(20+X2)X2 sinxi] - Kxi +X 2 ) 

2Md^ 


/ 


-[S + cosO^fi -;^f2) + 3cos;ifi 


+ 3 cos;r 2 ] + 3/|. (10c) 

There are two controlled parameters ix\,X 2 ) that are used to 
specify the shape of the swimmer. When the periodic motion 
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Fig. 4. (Color online) Time evolution of (a) the density profile and 

(b) the position of the COM of a three-body swimmer driven as in Eq. (11), 
with A = 7r/4 and y = oj = 0.025t“^ See the Supplemental Material for a 
movie of the dynamics in (a).^^ The vertical line in (b) indicates t = 350t. 

(c) Flow field J at t = 350t. The length of the vector is proportional to |7|, 
and the color represents Jx. The black arrows indicate the directions in which 
the joints are rotating, (d) Time evolution of the position of the COM of the 
three-body swimmer driven as in Eq. (13), with A = 7r/4, y = 2 x 10 “^t“^, 
to = 800t, and oj = 0.025t“^ 


of the swimmer is represented by a closed loop in (xi,X 2 ) 
space, the motion is neither time-reversal symmetric nor re¬ 
ciprocal,^^ and therefore, locomotion is not prohibited by the 
generalized scallop theorem. 

Figures 4(a) and 4(b) show the dynamics of the three-body 
swimmer for 

X\(,t) = -a|i - 2tanhyrsin^ yj, (11a) 

Xlif) = A tanh yf sin , (Hb) 


Since the trajectories of;^i(0 andA2(0 are circles in 
iXi^Xi) space, net locomotion is allowed. We can see from 
Figs. 4(a) and 4(b) that the three-body swimmer propels itself 
without exciting waves or vortices, in contrast to the two-body 
swimmer. The three-body swimmer travels by ^ 3^ for each 
stroke cycle. 

Intuitive understanding of self-propulsion is difficult, be¬ 
cause the flow field 

J = (12) 

Imi 

around the three-body swimmer is complicated, as shown in 
Fig. 4(c). Ai t = 350t, the position X of the COM of the 
swimmer is increasing [vertical line in Fig. 4(b)], and hence 
the superfiuid should have a momentum in the -v direction, 
due to conservation of the total momentum. The red (light 
gray) region in Fig. 4(c) has a large momentum in the -v di¬ 
rection, and this plays an important role in the locomotion of 
the swimmer. 

Figure 4(d) shows the position of the COM for swimmer’s 
shape given by 


Ai(0 

= A(l-2 c-^('-'»>%in2^), 

(13a) 

T2(0 

= sin yt, 

(13b) 


where A = 7r/4,y = 2 x10“^t“^, Iq = 800t, andcc) = 0.025 t“^. 
These functions begin at ai( 0) - ^/4 and a(2 (0) - 0, os¬ 
cillate while tjr ^ 800 + 300, and stop at ai(^) = ^/4 
and A 2 (<^) = 0- In this case, the swimmer travels only dur¬ 
ing the deformation, and the locomotion nearly stops after the 
shape stops changing. This behavior is also similar to the self¬ 
propulsion of a swimmer in a viscous Stokes fiuid. The small 
fluctuations of X and Y for f/r > 1200, shown in Fig. 4(d), are 
due to a non-adiabatic disturbance of the condensate. They 
vanish when the deformation of the swimmer is sufficiently 
slow. 

3. Swimming in multicomponent condensates 

We next consider the dynamics of a multicomponent EEC, 
in which there is a swimmer who also consists of BECs of 
different components. We will consider an immiscible three- 
component BEC. Suppose that a swimmer of components 1 
and 2 is surrounded by a sea of component 3. By controlling 
the scattering lengths by using the Eeshbach resonance, the 
size of the droplets of components 1 and 2 can be changed, 
that is, the shape of the swimmer can be controlled. 

The dynamics of the system can be studied by numerically 
solving the GP equation for a three-component BEC 

j 9 9 

^ ^ ( 14 ) 

k=l 


where A = 7r/4 and y = w = 0.025t“*. The factor tanhyr is where j = 1,2, and 3. The interaction coefficients gjt depend 
introduced to avoid infinite acceleration of the body at f = 0. on time in such a way that they always satisfy the immiscible 
For t » y"', Eq. (11) reduces X.ox\{t) - -coscot wAxiit) - 


4 





















J. Phys. Soc. Jpn. 


FULL PAPERS 


(a) 


t = 0 

t = 120x 

t = i50x r 

0 

0 

0 


' 30^ ' 

L 



Fig. 5. (Color online) Time evolution of (a) the density profile |j/r 3 p of 
component 3, and (b) the position X of the COM of the swimmer. See the 
Supplemental Material for a movie of the dynamics shown in (a).^^ 


condition: > gjjgkk, where j ^ k. They are changed as 


g 

gn = -(l-\-t2inhytcosojt). 

(15a) 

g 

g 22 = 2^^ tanhyr sin ojt), 

(15b) 


where y is introduced to avoid an initial disturbance of the 
system and g is a parameter that determines the units v^, 
and T. In the following calculation, we take y = 0.0 1t“^ and 
oj = 27rx0.002T“^ The other interaction coefficients are fixed 
as follows: gss = 1.5g and gn = g 23 = gi 3 = g- For t » 
the interaction coefficients in Eq. (15) become gn oc 1 + 
cos ojt and g 22 oc 1 + sin ojt, which do not have time-reversal 
symmetry. 

The initial state is the ground state with f = 

f \il/ 2 \^dr = lOO^o^^- The density 10 ^ 3 p of component 3 far 
from the center is a constant. The initial density profile 
of 103 p is shown in the leftmost panel of Fig. 5(a), where the 
left and right sides of the low-density region are occupied by 
components 1 and 2, respectively. By changing the interaction 
coefficients as in Eq. (15), the shape of the swimmer deforms, 
as shown in Fig. 5(a). We find from Fig. 5(b) that the swimmer 
obtains net locomotion in the -v direction without generating 
waves or quantized vortices. 

The locomotion of the swimmer in Fig. 5 may be under¬ 
stood qualitatively as follows. Suppose that the shapes of 
components 1 and 2 are approximated by circles with radii r\ 
centered at x\ and r 2 at V 2 , respectively, and that these circles 
touch each other, i.e., 

ri r 2 = X 2 - x\ > 0 . (16) 

When component 2 expands or shrinks, component 1 is 
pushed or pulled as vi ^ -r 2 , if the density ni of component 1 


is the same as the density of the surrounding fiuid of com¬ 
ponent 3. If ^1 > n^{n\ < ^ 3 ), the displacement of component 
1 is suppressed (enhanced), which may be approximated by 
xi ^ -r 2 + e{f - ri )f 2 to the first order of f - ri, where f is the 
radius for which m = and 6 > 0 is a constant. The same is 
also true when component 1 expands or shrinks, and then we 
have 


Xi -r 2 + ^(r - ri)r 2 - e{r - r 2 )h, (17a) 

■T 2 - h- e{r - r2)h + e{r - ri)r 2 , (17b) 

where the third terms on the right-hand side ensure Eq. (16). 
The motion of the COM of the swimmer thus becomes 

X = i(i:i -r ± 2 ) - ^(ri - h) - ^(r - r 2 )ri + e(r - ri)r 2 . (18) 

Assuming that ri - r cx cos ojt and r 2 - f oc sin ojt, the time- 

averaged COM in Eq. (18) becomes X oc -eoj, which agrees 
with the fact that the swimmer travels in the -v direction in 
Fig. 5. It has been confirmed from numerical simulations of 
Eq. (14) that the time-averaged velocity X is proportional to 
OJ, in agreement with the above result. 

The following facts have been confirmed numerically. In 
a manner similar to what is shown in Fig. 4(c), the self¬ 
propulsion ceases when the swimmer stops deformation. If 
the change in the interaction coefficients has time-reversal 
symmetry, no net locomotion is obtained, due to the gener¬ 
alized scallop theorem. For a larger oj, the swimmer emits 
excitations, in a manner similar to that shown in Fig. 3. How¬ 
ever, in this case, components 1 and 2 , that is, the body of the 
swimmer, are also emitted with waves and solitonic objects of 
component 3, and the swimmer gradually shrinks. 

4. Conclusions 

We have investigated the dynamics of deformable objects 
in superfiuids in order to answer the question: can we swim in 
a superfiuid? We began by examining the dynamics of artic¬ 
ulated bodies, as shown in Fig. 1. For a swimmer that con¬ 
sists of two elliptic bodies, no net locomotion is achieved 
when the shape of the swimmer changes sufficiently slowly 
(Fig. 2); this is consistent with the generalized scallop theo¬ 
rem. When the swimmer is deformed with a velocity compa¬ 
rable to or faster than the sound velocity, waves and solitonic 
excitations are emitted into the superfiuid, and the swimmer 
thus acquires momentum (Fig. 3). A swimmer consisting of 
three elliptic bodies can swim in a superfiuid without exciting 
waves (Fig. 4), since its motion can break time-reversal sym¬ 
metry (non-reciprocal^^). We next considered the dynamics 
of a three-component BEC, in which a swimmer comprising 
components 1 and 2 is surrounded by component 3. The shape 
of the swimmer can be changed by using the Feshbach reso¬ 
nance to control the interactions. We found that it is possible 
for such a swimmer to create self-propulsion without excit¬ 
ing waves; this is shown in Fig. 5. Thus, it has been shown 
that swimming in a superfiuid can be done smoothly or with 
splashing. 
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Although, for simplicity, we have considered only uniform 
2D systems, self-propulsion will also be possible in 3D super¬ 
fluids. However, spatial inhomogeneity, such as in an atomic 
gas EEC, may hinder self-propulsion. A solution may be 
possible in a system with ring geometry that has rotational 
symmetry instead of translational symmetry; in such a sys¬ 
tem, self-propulsion may be verifled by a swimmer traveling 
around the circumference of the ring. An interesting area for 
future research is to consider self-propulsion methods that are 
unique to quantum fluids, such as swimming by using quan¬ 
tized vortices or matter-wave interference. 
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